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ABSTRACT 

We  investigate  whether  ROSAT  observations  of  thermal  X-ray  emission  from  shell-type  supernova  rem¬ 
nants  (SNRs)  can  be  used  to  constrain  their  distances.  The  most  critical  assumption  is  that  the  initial  kinetic 
energy  (e0)  of  a  supernova  (SN)  explosion  is  a  constant  at  e0  =  1051  ergs;  however,  the  derived  distance  (D)  is 
only  weakly  dependent  on  e0  (D  oc  e£'4).  We  evaluate  this  technique  by  applying  it  to  SNRs  with  independent¬ 
ly  determined  distances,  and  our  initial  results  indicate  good  agreement.  We  conclude  that  the  ROSAT  all-sky 
survey  may  be  used  to  establish,  for  the  first  time,  a  set  of  good  distance  estimates  to  a  large  number  of 
extended,  shell-type  SNRs.  The  energy  range  of  ROSAT  is  well  suited  for  this  purpose,  since  most  shell-type 
SNRs  have  thermal  X-ray  spectra  which  peak  within  the  ROSAT  PSPC  energy  range  of  0. 1-2.4  keV. 

Subject  headings:  supernova  remnants  —  X-rays:  general 


1.  INTRODUCTION 

SNRs  are  important  in  Galactic  astrophysics  because  they 
are  the  major  source  of  energy  input  to  the  interstellar  medium, 
serve  as  the  site  of  cosmic-ray  acceleration  and  play  a  role  in 
triggering  star  formation.  Unlike  H  n  regions,  whose  distances 
can  be  constrained  kinematically,  Galactic  SNRs  have  poorly 
determined  distances.  Of  the  currently  known  ~  ISO  Galactic 
remnants,  only  ~  15%  have  distance  determinations  and  only 
a  limited  subset  of  these  are  likely  to  be  reliable  (Green  1984). 
Most  methods  of  establishing  distances  to  SNRs  are  either 
difficult  and  uncertain  (H  I  absorption)  or  are  too  inexact  to  be 
useful  (e.g.,  the  “  E-D  ”  relation)  (Green  1984).  We  describe  how 
soft  X-ray  observations,  such  as  those  made  by  the  German 
X-ray  observatory  ROSAT  (Triimper  1983),  can  be  used  to 
estimate  distances  to  many  shell-type  SNRs.  We  have  already 
applied  this  procedure  to  G326.3— 1.8  (Kassim,  Hertz,  & 
Weiler  1993),  and  we  apply  it  here  to  a  number  of  SNRs  with 
independently  established  distances  in  order  to  test  the  pro¬ 
cedure’s  validity. 

2.  THEORY 

The  Sedov  equations,  which  describe  the  dynamics  of  shell- 
type  SNR  expansion  during  the  adiabatic  phase,  can  be  com¬ 
bined  (Kassim  et  al.  1993)  to  solve  for  a  “  Sedov  ”  distance  Ds  to 
the  SNR  in  kiloparsecs  as 

Ds  =  8.7  x  106eg-4P(A£,  t)°-20-°-6F;oo-2T-0-4  .  (1) 

Here  e0  is  the  initial  energy  of  the  SN  explosion  (thermal  plus 
kinetic)  in  units  of  10s  1  ergs,  0  is  the  observed  angular  diameter 
of  the  SNR  shell  in  arcminutes,  Fx0  is  the  measured  X-ray  flux 
corrected  for  interstellar  (and/or  intergalactic)  absorption  in 
units  of  ergs  s-1  cm-2,  and  T  is  the  measured  thermal  tem¬ 
perature  of  the  X-ray-emitting  gas  in  K.  The  tabulated  analyti¬ 
cal  function  P(AE,  T)  in  units  of  ergs  cm3  s~ 1  describes  the 
power  emitted  by  hot  electrons  in  a  low-density  plasma  via 
free-free  emission  (Tucker  &  Koren  1971a,  b)  and  is  a  function 
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of  both  the  energy  band  of  the  emission  and  the  temperature  of 
the  plasma.  The  uncertainties  in  determining  Ds  are  dominated 
by  the  assumptions  inherent  in  applying  equation  (1)  and  not 
by  measurement  errors  in  the  temperature,  interstellar  absorp¬ 
tion,  and  flux  as  determined  from  ROSAT  PSPC  X-ray 
spectra,  as  was  often  the  case  with  pre-R  OS/IT  X-ray  observa¬ 
tions. 

The  key  assumptions  of  the  Sedov  analysis  reflected  by 
equation  (1)  are  (1)  the  SNR  shell  is  in  the  adiabatic  expansion 
phase;  (2)  the  measured  X-ray  temperature  gives  a  reliable 
estimate  of  the  SN  shock  velocity;  and  (3)  e0  =  1  (in  units  of 
10s  1  ergs)  is  approximately  correct  for  all  SNe.  Assumption  (1) 
is  the  only  one  which  can  be  checked  self-consistently.  A  calcu¬ 
lation  of  the  swept-up  mass  (MSJ  can  be  used  to  test  whether  a 
SNR  has  passed  the  free-expansion  phase  and  entered  the  adia¬ 
batic  expansion  phase  (Kassim  et  al.  1993),  and  relations  have 
been  developed  (Cox  1972;  McKee  1982;  Cioffi,  McKee,  & 
Bertschinger  1988)  to  calculate  if  the  SNR  has  entered  either 
the  pressure-driven  or  radiative  phases. 

Unfortunately,  neither  assumptions  2  nor  3  can  be  firmly 
established  on  either  theoretical  or  observational  grounds. 
Calculations  of  Type  I  SNe  show  that  e0  only  varies  over  the 
range  0.9-1. 5  (Nomoto,  Thielemann,  &  Yokoi  1984),  but  there 
are  no  generally  accepted  ways  of  estimating  the  energy  of 
Type  II  events.  With  the  total  energy  released  in  the  gravita¬ 
tional  collapse  of  a  star  during  a  SN  explosion  of  ~  1053  ergs, 
the  assumed  thermal  plus  kinetic  energy  e0  of  1051  ergs  is  only 
~  1%  of  the  total.  Since  the  efficiency  of  transforming  gravita¬ 
tional  energy  into  thermal  and  kinetic  forms  is  not  well  under¬ 
stood  and  could  be  highly  variable,  assumption  3  may  be  only 
roughly  correct.  Assumption  2  is  related  to  complex  issues 
involving  the  degree  of  thermalization  of  the  kinetic  energy  in 
the  shock,  the  extent  of  the  equilibrium  between  ion  and  elec¬ 
tron  energies,  and  the  equilibration  of  the  electron  velocity 
distribution  and  ionization  nonequilibrium.  These  areas  are 
directly  related  to  poorly  determined  theories  of  SNe  and 
strong  shocks.  Also,  the  detection  of  X-rays  from  cooler  gas 
behind  the  shock  front  and  in  the  remnant  interior  may  lead  to 
a  measured  X-ray  temperature  which  is  less  than  the  shock 
temperature  (Rappaport  et  al.  1974),  and  hence  an  underesti¬ 
mate  of  the  SN  shock  velocity.  We  also  note  that  the  assump¬ 
tion  of  solar  abundances  in  the  line-dominated  X-ray  spectra 
could  lead  to  errors  in  the  X-ray  temperature  when  measured 
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with  a  low  spectral  resolution  device  like  the  ROSAT  PS  PC 
(Hughes  1994). 

3.  APPLICATION 

Fortunately,  we  can  conduct  an  empirical  test  of  the  validity 
of  the  “  Sedov”  distance  as  expressed  in  equation  (1)  by  apply¬ 
ing  it  to  ROSAT  data  for  shell- type  SNRs  with  relatively  well- 
established  independent  distance  estimates.  Table  1 
summarizes  the  data  for  our  sample.  All  values  of  T  and  Fx0 
(col.  [2]-[3])  are  from  ROSAT  observations  (energy  range 
~0.1-2.4  keV)  as  referenced  in  column  (9),  except  for  the  LMC 
SNR  N49  (Gordon  et  al.  1993).  The  mainly  radio-determined 
angular  sizes  (0;  col.  [4])  are  obtained  from  Green’s  SNR 
catalog  (Green  1990, 1991,  1994)  for  the  Galactic  SNRs,  other¬ 
wise  as  referenced  in  column  (9).  The  radio  morphology  is 
known  to  be  a  good  tracer  of  the  blast  wave,  while  the  X-ray 
morphology  can  be  centrally  condensed  even  for  shell-type 
SNRs  (Dickel  et  al.  1994;  Jones,  Smith,  &  Angellini  1993;  Rho 
et  al.  1994).  The  swept-up  mass  (Msu,  col.  [5])  is  determined  by 
solving  the  Sedov  equation  (see,  e.g.,  Kassim  et  al.  1993) 

Msu  =  3.6  x  10~9P(AE,  T)-°  s0i  :>F%sD2-5  ,  (2) 

with  the  assumption  that  the  distance,  D,  is  given  by  the  inde¬ 
pendent  distance  estimate  (Dh  col  [7]);  no  assumption  about 
the  explosion  energy  is  required. 

The  Sedov  distances  ( Ds ,  col.  [6])  are  obtained  from  equa¬ 
tion  (1).  The  independent  distances  (Dt,  col.  [7])  are  obtained 
from  the  Green  SNR  catalog  (Green  1990,  1991,  1994)  when 
available,  otherwise  as  referenced  in  column  (9).  The  ratio 
Ds/Dt  (col.  [8])  is  also  listed,  and  the  two  distance  estimates  are 
plotted  against  one  another  in  Figure  1  for  all  of  the  remnants 
except  G18.9— 1.1,  which  is  probably  not  in  the  Sedov  phase 
(see  below).  The  error  bars  represent  statistical  uncertainties 
only  and  are  based  on  propagating  measurement  errors  for  T, 
0,  Fx0,  and  D(  where  available;  otherwise  we  use  canonical 
uncertainties  of  10%  for  0,  25%  for  T  and  Du  and  the  smaller 
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INDEPENDENT  DISTANCE  (kpc) 

Fig.  1. — Plot  of  the  X-ray  determined  Sedov  distance  Ds  vs.  the  indepen¬ 
dently  determined  distance  Dt  for  seven  Galactic  and  three  extragalactic  shell- 
type  SNRs.  Error  bars  represent  statistical  uncertainties  only  (see  text).  The 
data  are  clearly  well  fitted  by  the  dotted  line  representing  Ds-  Dt. 

of  20%  and  2  x  10  12  ergs  cm  “ 2  s  “ 1  for  Fx0.  Note  that  the 
measurement  errors  from  ROSA  T -determined  Sedov  distances 
are  smaller  than  the  typical  uncertainty  in  H  i  distance  deter¬ 
minations  for  Galactic  SNRs.  Figure  1  shows  clearly  that 
Ds  =  Dt  is  consistent  to  within  the  errors. 

4.  DISCUSSION 

With  the  single  exception  of  G18.9— 1.1,  discussed  below,  all 
SNRs  listed  in  Table  1  have  Msu  >  20  M0,  implying  that  they 


TABLE  1 

Observed  and  Derived  Radio  and  X-Ray  Parameters  tor  1 1  Supernova  Remnants 


kT“  iV  e*  Afsu'  Dsd  D t' 

SNR  (keV)  (ergs  s  1  cm  2)  (')  (Me)  (kpc)  (kpc)  DJD,  Reference 

(0  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9) 


G1 8.9 -1.1  .  0.43  2.8E-11  33  4.4  8.4  2.0  4.2  1 

G160.9  +  2.6  .  0.61  1.0E-10  130  64  2.5  2.0f  1.3  2 

G156.2+5.7  .  0.50  1.9E-10  108  115  2.7  2.5  1.1  3 

G34.7— 0.4 .  2.8  2.5E-10  31  46  2.3  3.tf  0.8  4 

G132.7+ 1.3  .  0.16  9.3E-10  80  191  4.2  3.0f  1.4  2 

G296. 1—0,7  .  0.20  6.3E-10  33  98  6.6  4.0  1.7  5 

G53.6— 2.2 .  0.27  1.9E-10  28  174  7.8  6  7f  12  2 

G33.6+0.1  .  1.3  4.3E-11  10  62  9,4  10f  0.9  6 

N49 .  0.65  8.0E-11  1.1  119  48  50s  1.0  7 

0540-69.3  .  5.5  2.3E-12  0.9  23  39  50e  0.8  8 

013022+30233  .  0.33  4.8E-13  0.1  272  758  840b  0.9  7 


a  All  ROSAT  determined  except  for  the  LMC  SNR  N49  (Gordon  et  al.  1993). 

h  Radio-determined  (Green  1990, 1991, 1994)  for  Galactic  SNRs,  otherwise  as  referenced  in  col.  (9). 

c  All  swept-up  masses  are  determined  from  the  Sedov  equations  assuming  £>,  is  the  correct  distance  (see  text). 

d  Sedov  distance  from  eq.  (1). 

e  As  noted;  otherwise  as  referenced  in  (9). 

f  From  the  Green  SNR  catalog  (Green  1990, 1991, 1994). 

*  Adopted  distance  to  the  LMC  (Panagia  et  al.  1991). 
h  Adopted  distance  to  M33  (Freedman,  Wilson,  &  Madore  1991). 

References.— (1)  Aschenbach  et  al.  1991;  (2)  Saken  et  al.  1994;  (3)  Pfelfermann,  Aschenbach,  &  Predehl  1991;  (4) 
ROSAT  Archive;  (5)  Hwang  &  Markert  1994;  (6)  Seward  &  Velusamy  1993;  (7)  Gordon  et  al.  1993;  (8)  Seward  & 
Harnden  1994. 
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have  passed  the  free-expansion  phase.  Also,  we  have  calculated 
the  radii  at  which  the  SNRs  in  Table  1  are  expected  to  enter  the 
pressure  driven  and  radiative  phases  and,  under  the  criteria  of 
Burrows  et  al.  (1993)  with  e0  =  1,  none  has  yet  evolved  beyond 
the  adiabatic  phase. 

While  the  correlation  between  Ds  and  D,  appears  well 
defined  for  all  SNRs  in  Figure  1,  Table  1  includes  one  object, 
G  18.9  — 1.1,  which  illustrates  that  the  technique  is  not  applic¬ 
able  to  SNRs  which  have  not  fully  evolved  to  the  Sedov  stage. 
An  independent  distance  estimate  of  ~  2  kpc  from  H  i  observa¬ 
tions  has  been  used  by  Aschenbach  et  al.  (1991)  with  the  Sedov 
assumption  to  calculate  e0  =  0.05.  However,  the  estimated 
swept-up  mass  for  this  SNR  is  only  4.4  M0.  Since  existing 
models  for  SN  ejected  masses  range  from  ~  1-5  M0  (Woosley 
et  al.  1994),  it  is  quite  likely  that  G18.9  — 1.1  is  still  in  transition 
from  the  free-expansion  to  the  Sedov  stage,  and  therefore  the 
Aschenbach  et  al.  (1991)  calculation  of  the  explosion  energy 
may  not  be  valid.  For  the  same  reason,  our  distance  estimating 
technique  probably  does  not  apply  to  G  18.9  — 1.1,  and  we  have 
excluded  it  from  Figure  1. 

The  average  value  of  Ds/D,  for  the  remaining  10  SNRs 
included  in  Figure  1  is  1.1,  with  a  standard  deviation  of  ~0.3. 
To  within  the  errors,  this  implies  that,  for  SNRs  in  the  Sedov 
phase,  our  X-ray-based  distance  estimates  are  in  agreement 
with  distances  obtained  by  other  methods  and  the  technique 
can  likely  be  applied  to  the  many  SNRs  observed  with  ROSAT 
for  which  independent  distance  estimates  are  not  available. 


5.  SUMMARY 

Our  initial  test  comparing  ROSAT  derived  “Sedov”  dis¬ 
tances  with  independent  distance  estimates  for  shell-type  SNRs 
indicates  the  modern  X-ray  measurements  are  a  powerful  tool 
for  obtaining  distances  to  large  numbers  of  shell-type  SNRs. 
While  several  key  assumptions  cannot  be  rigorously  justified 
on  theoretical  grounds  at  this  time,  the  empirical  conclusion  is 
that  the  method  offers  a  useful  distance  estimate  in  the  absence 
of  other  data.  High-resolution  X-ray  spectra  of  SNRs,  such  as 
those  obtainable  with  the  Japanese  ASCA  ( Astro  D)  satellite, 
will  test  our  ability  to  measure  plasma  temperatures  in  line 
dominated  plasmas  with  a  low-resolution  spectrometer,  such 
as  the  ROSAT  PSPC.  As  the  X-ray  parameters  become  avail¬ 
able  for  many  more  SNRs  from  the  ROSAT  all-sky  survey, 
further  tests  of  the  method  will  better  define  its  accuracy  and 
utility. 
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